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Preface 

This  analysis  was  conducted  at  the  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  CE,  Hydraulics  Laboratory  (HL) ,  from  September 
1976  to  December  1978.  This  effort  was  initiated  under  work  unit 
CWIS  31402,  '’Methods  of  Enhancing  Water  Quality,”  of  the  Reservoir  Water 
Quality  Research  Program  sponsored  by  the  Office,  Chief  of  Engineers. 

The  effort  was  completed  and  the  report  published  within  the  Environmen¬ 
tal  and  Water  Quality  Operational  Studies  (EWQOS)  under  EWQOS  Work 
Unit  31604  (III  A.l). 

Dr.  J.  L.  Mahloch  was  Program  Manager  for  EWQOS-  Messrs*  H.  B. 
Simmons,  Chief  of  the  HL,  and  J.  L.  Grace,  Jr.,  Chief  of  the  Hydraulic 
Structures  Division,  directed  the  effort.  Mr.  Steven  C.  Wilhelms  con¬ 
ducted  the  study  and  prepared  the  text  under  the  direct  supervision  of 
Mr.  D.  G.  Fontane,  former  Chief  of  the  Reservoir  Water  Quality  Branch 
(Physical).  Dr.  D.  R.  Smith,  Chief  of  the  Reservoir  Water  Quality  Branch, 
assisted  in  the  preparation  and  review  of  this  report.  Mr.  M.  E.  Neumann 
assisted  in  the  mathematical  analysis  of  the  data.  Data  were  furnished 
by  Mr.  Glenn  Drummond  and  other  personnel  of  the  Ohio  River  Division  and 
by  personnel  from  the  Vicksburg  District,  CE. 

COL  John  L.  Cannon,  CE,  and  COL  Nelson  P.  Conover,  CE,  were 
Commanders  and  Directors  of  WES  during  this  effort.  Mr.  Fred  R.  Brown 
was  Technical  Director. 

This  report  should  be  cited  as  follows: 

Wilhelms,  Steven  C.,  and  Smith,  Dennis  R.  1981. 

"Reaeration  Through  Gated-Conduit  Outlet  Works," 

Technical  Report  E-81-5,  U.  S.  Army  Engineer  Water¬ 
ways  Experiment  Station,  CE,  Vicksburg,  Miss. 
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Conversion  Factors,  U.  S.  Customary  To  Metric  (SI) 
Units  of  Measurement 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  (SI)  units  as  follows: 


_ Multiply _ 

cubic  feet  per  second 
Fahrenheit  degrees 
feet 

foot-pounds  ( force) 
pounds  (force) 
square  feet 


By_ 

0.02832 

5/9 

0. 3048 
1.355818 
4.44822 
0.09290304 


_ To  Obtain _ 

cubic  metres  per  second 
Celsius  degrees  or  Kelvins* 
metres 

metre-newtons  or  joules 

newtons 

square  metres 


*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  read¬ 
ings,  use  the  following  formula:  C  ■  (5/9) (F  -  32).  To  obtain 
Kelvin  (K)  readings,  use:  K  -  (5/9) (F  -  32)  +  273.15. 
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REAERATION  THROUGH  GATED-CONDUIT  OUTLET  WORKS 


Introduction 


1.  It  is  often  possible  to  enhance  the  water  quality  of  releases 
from  reservoirs  by  employing  outlet  works  that  aerate  the  flow.  This 

is  especially  important  if  the  water  entering  the  outlet  works  is  low  in 
dissolved  oxygen  (D.O.),  since  many  kinds  of  aquatic  life  cannot  survive 
in  a  low  D.O.  environment.  To  design  structures  that  effectively 
aerate  the  relase  water,  predictive  techniques  are  required  that  quantify 
the  gas  transfer  that  results  with  flow  through  a  particular  structure 
of  interest.  Few  predictive  techniques  have  been  developed  or  applied  to 
hydraulic  structures. 

2.  The  development  of  generalized  equations  applicable  to  arbi¬ 
trary  structural  designs  is  extremely  difficult.  Reaeration  during  flow 
through  hydraulic  structures  depends  upon  molecular  diffusion  and  the 
magnitude  of  dispersion  due  to  turbulent  mixing.  In  most  instances,  tur¬ 
bulence  is  the  physical  mechanism  that  dominates  the  rate  of  gas  transfer. 
Numerical  modeling  of  turbulence  to  accurately  describe  gas-transfer 
processes  in  self-aerated  flows  requires  a  data  base  and  modeling 
techniques  in  excess  of  those  currently  available.  Reaeration  occurring 
during  flow  through  hydraulic  structures  is  currently  predicted  on  a 
very  limited  basis  by  use  of  semiempirical  or  empirical  mathematical 
models  of  the  reaeration  process.  For  any  of  these  models,  the  coeffi¬ 
cients  must  be  determined  experimentally  and  the  models  compared  to 
determine  the  most  applicable.  Above  all,  the  models  should  not  be 
applied  indiscriminately  outside  of  the  flow  regime  for  which  the  coeffi¬ 
cients  were  developed.  It  is  imperative  to  investigate  the  adequacy  of 
the  respective  models  for  predicting  reaeration  in  hydraulic  structures. 

3.  Two  models  were  analyzed  and  evaluated  regarding  reaeration 
occurring  through  outlet  works:  (a)  Energy  Dissipation  Model  (EDM) 
(Tsivoglou  and  Wallace  1972)  and  (b)  Deficit  Ratio  Model  (DRM)  (Holler 
1970),  These  two  were  chosen  because  they  are  relatively  simple  mathe¬ 
matical  models  that  predict  the  gas  transfer  as  a  function  of  the  total 
energy  dissipated. 
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4.  The  objectives  of  this  study  were  as  follows: 

a.  To  analyze  and  document  existing  D.O.  data  for  gated- 
conduit  release  structures. 

To  provide  a  qualitative  and  quantitative  evaluation  of 
the  gas-transfer  characteristics  of  existing  structures 
of  this  type. 

£,  To  test  the  applicability  of  the  EDM  and  DRM  for  predict¬ 
ing  oxygen  uptake  through  reservoir  outlet  works. 


Scope 

5.  This  report  addresses  reservoir  outlet  works  consisting  of 
intake  structures,  conduits,  and  stilling  basins  (Figure  1).  Analysis 
of  D.O.  data  from  several  U.  S.  Army  Corps  of  Engineers  projects  is  re¬ 
ported  herein.  Drawings  of  the  projects  are  presented  in  Plates  1-12. 

All  of  these  outlet  works  exhibited  similar  hydraulic  conditions  inas¬ 
much  as  there  was  free-surface  flow  through  the  conduit  and  the  outlet 
portal  was  not  submerged  by  the  tailwater. 

Approach 

6.  The  data  available  for  analysis  were  observed  D.O.  and  tempera¬ 
ture  profiles  in  the  lake  and  downstream  D.O.  and  temperature  (Plates 
13-38).  Prediction  of  the  D.O.  and  temperature  entering  the  structure 
was  required.  A  generalized  selective  withdrawal  technique  (Bohan  and 
Grace  1973),  "SELECT,"  developed  at  the  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  CE,  was  verified  for  applicability  to  the 
hydraulic  structures  by  predicting  release  temperatures  and  comparing 
the  predictions  with  the  observed  release  temperatures.  This  was  possi¬ 
ble  because  there  is  negligible  heat  transfer  between  the  release  water 
and  structure.  SELECT  was  then  used  to  predict  D.O.  content  of  the 
water  entering  the  hydraulic  structure.  These  predictions  were  used  in 

a  qualitative  evaluation  of  the  reaeration  occurring  through  the 


structures  and  as  input  to  the  EDM  and  DRM  (Appendix  A). 
Energy  Dissipation  Model 


7.  The  E DM ,  which  was  developed  for  streams,  states  that  the  re¬ 
aeration  rate  coefficient  is  proportional  to  the  rate  of  energy  dissipa¬ 
tion,  The  model  is  expressed  mathematically  as 


where 


(1) 


=  reaeration  or  gas-transfer  coefficient,  base  e  ,  per  sec 

c  =  escape  coefficient,  per  ft*  of  energy  loss 

AE  =  energy  expended  from  upstream  point  to  downstream  point, 
ft-lb  per  lb  of  water 

t^  =  time  of  flow  from  upstream  to  downstream,  sec 

The  total  energy  dissipated  between  upstream  and  downstream  points  is 

the  sum  of  the  irreversible  losses  in  kinetic  and  potential  energy, 

which  can  be  expressed  as 

2  2 
V  -  V 

AE  =  ~~ — -  +  (Hx  +  ZL)  -  (H2  +  Z2)  (2) 


where  the  variables  are  defined  by  (Tsivoglou  and  Wallace  1972) 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (SI)  units  is  presented  on  page  3. 


where 


,  ^2  =  vel°c!^  respective  locations,  ft/sec 
Hi*  =  depth  of  flow  at  respective  locations,  ft 

Z  ,  Z  =  elevation  of  respective  locations  above  some  datum,  ft 

1  1  2 
g  =  acceleration  due  to  gravity,  ft/sec 

Ah  =  water  surface  elevation  change  from  upstream  to 
downstream,  ft 

In  the  hydraulic  structures  analyzed,  the  velocity  in  the  lake  was 

approximately  zero,  and  the  velocity  in  the  exit  channel  was  small 

relative  to  its  contribution  to  total  head.  Thus,  the  term 


2g 


is  negligible  when  compared  to  the  elevation  difference.  Therefore, 

AE  can  be  closely  approximated  by  Ah  the  water  surface  elevation 
change  from  upstream  to  downstream  expressed  in  feet: 

AE  ~  Ah  =  (H2  +  Z2)  -  (H2  +  Z2)  (3) 

8.  Assuming  no  chemical  or  biological  oxygen  consumption  during 
flow  through  the  structure,  the  downstream  D.O.  concentration  can  be 
estimated  from  a  first  order  reaction  equation  (Streeter  and  Phelps 
1925) 


Df  =  Di  exp  (-K2tf) 


(4) 


where  equals  D.O.  deficit  downstream  and  upstream,  re¬ 

spectively,  mg/ it.  The  oxygen  deficit  is  defined  as 


where 


D 


sat 


C 

ac 


t 


D  = 


sat 


C 


act 


D.O.  deficit,  mg /l 

saturation  concentration  of  D.O.  ,  mg/ Jt 
actual  concentration  of  D.O.  ,  mg /£, 
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Combining  Equations  1,  3,  and  4  yields 


exp  (-cAh) 


(5) 


the  predictive  EDM. 

9.  Tests  on  streams  and  rivers  (Tsivoglou  and  Neal  1976)  using 
a  gaseous  tracer  (Tsivoglou  et  al.  1968)  indicated  that  the  escape 
coefficient  for  a  particular  reach  was  dependent  on  the  streamflow  rate. 
The  stream  observations  were  divided  into  groups  according  to  discharge. 
The  first  group  had  discharges  less  than  10  cfs.  The  second  group 
ranged  in  discharges  up  to  500-750  cfs.  A  third  group  was  considered 
(>750  cfs),  but  insufficient  data  were  available  for  analysis.  This 
information  implies  that  a  single  escape  coefficient  might  be  used  to 
predict  the  D.O.  uptake  occurring  in  this  type  of  structure  by  similar 
groupings  of  flow  rates.  Most  flows  in  these  projects  were  within  the 
range  of  the  second  group,  although  a  few  were  in  the  third  group. 

None  were  in  the  first  group. 

30.  The  escape  coefficient  for  gated-conduit  structures  was  un¬ 
known,  thus  gas-transfer  tests  were  conducted  at  the  outlet  works  of 
Enid  Lake,  Mississippi,*  using  a  gaseous  tracer  technique  (Tsivoglou 
et  al.  1968).  These  tests  indicated  a  value  of  0.045  per  ft  as  an  over¬ 
all  escape  coefficient  for  the  Enid  structure.  This  estimate  of  the 
escape  coefficient,  after  adjustment  for  the  ambient  water  temperatures 
by 


=  c 


20 


1.022 


(T-20) 


where  c^  ,  c0g  are  escape  coefficients  at  T°C  and  20°C,  respectively, 
was  used  in  the  predictive  EDM  (Equation  5)  to  estimate  oxygen  uptake 
occurring  in  several  hydraulic  structures  consisting  of  an  intake  struc¬ 
ture,  conduit,  and  stilling  basin. 


*  Unpublished  data,  C.  H.  Tate,  Jr.,  August  1978,  U.  S.  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  Mississippi. 
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Deficit  Ratio  Model 

11.  Using  the  "two  film"  theory  of  mass  transfer,  the  DRM  can  be 
expressed  (Holler  1970) 


where 

r  =  deficit  ratio 

=  liquid  film  coefficient 
a  =  air/liquid  interface  area 
t  =  time  of  contact 
v  =  volume  of  liquid 

12.  In  a  channel  the  average  transit  time  can  be  approximated  in 
terms  of  the  channel  volume  and  the  average  flow  rate: 


where 

Q  =  flow  rate 
t  =  time  of  flow 
Substituting  in  Equation  6, 


Using  a  Maclaurin!s  series  to  expand  Equation  7  gives 


n=l 


13.  It  has  been  theorized  (Holler  1970)  that  the  air  liquid 
interface  area  a  can  be  expressed  as  a  function  of  momentum  change 
occurring  in  a  hydraulic  jump  or  when  a  jet  impinges  on  a  surface.  By 
substituting  a  momentum  function  for  area  and  dropping  higher  order 
terms,  Equation  8  can  be  expressed  as 
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r  55  1  +  f  ( AV) 


where  f(AV)  is  some  function  of  the  velocity  change.  In  Holler's  (1970) 

2 

work,  the  deficit  ratio  was  experimentally  related  to  AV  and  yielded 
the  DRM  predictive  model 


r  =  1  +  B(AV)2 


(9) 


where 

AV  =  change  in  velocity,  ft /sec 

2  2 

6  =  coefficient,  sec  /ft 

For  flow  over  a  sharp-crested  weir,  kinematics  can  be  used  to  relate  the 
velocity  of  the  jet  to  fall  height  by  converting  potential  energy  to 
kinetic  energy  KE  of  the  jet. 


When  the  jet  impinges  on  a 
in  a  velocity  change.  The 


flat  plate,  the  momentum  change  is  reflected 
velocity  change  along  the  axis  of  the  jet  is 


AV  =  V 


since  there  is  no  velocity  along  the  direction  of  the  jet  after  impact. 
The  potential  energy  is  converted  to  kinetic  energy  and  a  large  part  of 

the  KE  is  dissipated  in  the  impact.  Holler  (1970)  substituted  h 
2 

for  AV  for  the  experimental  work  performed  on  sharp-crested  weirs. 

In  hydraulic  structures,  the  change  of  potential  energy  to  kinetic 
energy  and  subsequent  dissipation  is  analogous,  even  though  the  energy 
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is  not  dissipated  in  a  single  impact.  Thus,  Ah  ,  the  difference  in  pool 

2 

and  tailwater  elevations,  was  used  as  an  estimate  for  AV  in  Equa¬ 
tion  9  with  the  constant  2g  taken  into  account  in  the  coefficient  . 


Holler  determined  experimentally  that  8^  was  0.065  per  ft  for  dis¬ 
charges  over  low-head  submerged  tainter  gate  spillway  structures. 
Therefore , 


r  -  1  +  0.065Ah  (10) 

Equation  10  was  used  to  predict  D.O.  uptake  through  the  same  structures 
used  in  the  EDM  analysis.  Observed  D.O.  data  at  the  structures  were 
compared  to  predicted  D.O. 


Results 


14.  Plates  13-38  document  the  data  used  in  this  analysis  of 
reaeration  through  selected  outlet  works.  D.O.  and  temperature  profiles 
for  the  lakes  are  plotted,  and  observed  downstream  D.O.  and  temperature 
measurements  are  given.  The  Beltzville  data  are  documented  in  Hart  and 
Wilhelms  (1977).  Other  hydraulic  and  predicted  data  are  presented  in 
Appendix  A.  Release  water  temperatures  were  predicted  with  SELECT. 
Figure  2  shows  the  SELECT-predicted  release  temperature  plotted  against 
the  observed  release  temperature.  The  standard  deviation  (Miller  and 
Freund  1977)  for  SELECT-predicted  temperatures  was  +2.0°C. 

15.  The  data  from  all  the  projects  indicated  that  reaeration 
through  the  structures  was  sufficient  to  raise  the  downstream  D.O.  to 
approximately  80-100  percent  of  saturation  regardless  of  the  flow  and 
initial  D.O.  content.  Figure  3  shows  change  of  percent  saturation 
plotted  against  initial  saturation.  The  upstream  deficit  D^  was 
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Figure  2.  Predicted  temperature  versus  observed  temperature 
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Figure  3.  Change  in  D.O.  saturation  versus  initial  saturation 


predicted  with  SELECT.  Final  saturation  was  computed  from  observed 
data. 

16.  The  EDM  predictive  model  was  used  to  estimate  the  downstream 
deficit  .  Predictions  made  with  this  model  were  compared  to  ob¬ 

served  data  (Figure  4).  The  predicted  data  varied  with  a  standard  error 
(Miller  and  Freund  1977)  of  +0.7  mg/A.  D.O.  predictions  were  also  made 
with  the  DRM  and  were  compared  to  observed  downstream  D.O.  These 
comparisons  are  presented  in  Figure  5.  Standard  error  for  these  pre¬ 
dictions  was  +0.9  mg/fc. 


Conclusions 


17.  All  the  projects  analyzed  indicated  that  releases  reaerated 

to  approximately  90  percent  of  saturation.  A  common  minimum  release  D.O. 
objective  as  established  by  some  state  standards  is  5.0  mg/£.  Data 
clearly  indicated  that  these  gated-conduit  structures  successfully  met 
this  standard. 

18.  Most  D.O.  predictions  made  with  the  EDM  using  an  escape 
coefficient  c^q  *  0.045  per  ft  compared  favorably  with  observed  D.O. 

The  agreement  between  observed  and  predicted  data  was  very  good,  con¬ 
sidering  the  inherent  variability  in  the  data  and  the  hydraulic  differ¬ 
ences  among  the  structures.  It  is  probable  that  the  escape  coefficients 
vary  slightly  from  project  to  project  and  with  flow  for  a  particular 
structure  because  of  differences  in  levels  of  turbulence,  geometry,  or 
other  factors  that  affect  the  hydraulics  and  hydrodynamics  of  the  flow 
through  the  structure. 

19.  Osing  a  C£q  of  0.045  per  ft  for  gated-conduit  structures 
with  the  range  of  flow  rates  investigated  is  consistent  with  river  and 
stream  investigations  (Tsivoglou  and  Neal  1976) ,  which  have  indicated 
an  average  c^q  of  0.054  per  ft  for  flow  rates  similar  to  those  en¬ 
countered  in  these  projects.  It  was  concluded  from  those  studies  that 
as  the  flow  increases,  the  c^q  should  be  adjusted  downward  toward  a 
Uniting  value  of  approximately  0.025  to  0.030  per  ft. 

20.  D.O.  predictions  made  with  the  DRM  (Figure  5)  tend  to  be 
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Figure  5.  Holler's  (1970)  DRM-predi  cted  D.O.  versus  observed  D.O. 


lower  than  observed  D.O.  levels.  This  implies  that  the  coefficient  8' 
in  the  DRM,  which  was  developed  for  low  heads,  should  be  larger  for 
high-head  structures  to  increase  predicted  downstream  D.O.  concentration 
The  analysis  of  data  obtained  from  the  tracer  studies  at  Enid  Lake*  re¬ 
sulted  in  a  coefficient  of  0.200  per  ft.  Predictions  made  with  the 

DRM  with  this  coefficient  are  shown  in  Figure  6.  The  standard  error 
(Miller  and  Freund  1977)  of  prediction  was  +0.8  mg/fl,  as  compared  to 
0.9  mg/ l  obtained  with  8^  of  0.065  per  ft  developed  for  low-head 
structures . 

21.  The  above  results  indicate  that  there  is  very  little  differ¬ 
ence  in  the  standard  errors  for  the  EDM  model  (calibrated  with  the  Enid 
outlets  work  data),  the  DRM  calibrated  with  low-head  data,  and  the  DRM 
calibrated  with  the  Enid  data.  This  implication  is  somewhat  misleading, 
since  there  are  fundamental  differences  in  the  models.  The  differences 
can  be  clearly  delineated  by  comparing  plots  of  the  respective  equations 

22.  In  Figure  7  the  DRM  developed  by  Holler  (1970)  for  low-head 
submergible  tainter  gates  can  be  compared  with  the  EDM  calibrated  with 
the  Enid  data.  The  DRM  has  been  inverted  to  provide  consistency  in  the 
comparison.  The  models  predict  essentially  identical  deficit  ratios  for 
heads  less  than  approximately  20  ft  and  consequently  are  equally  valid. 
However,  above  20  ft,  the  differences  in  predicted  deficit  ratios  in¬ 
crease.  The  difference  at  larger  heads  is  inherent  in  the  derivation 

of  the  respective  models.  Both  assume  a  first  order  process  that  is 
normally  expressed  mathematically  as  an  exponential  process.  However, 
the  DRM  employs  a  truncated  Maclaurin's  series  expansion  to  describe  the 
exponential  process.  As  a  result,  the  DRM  is  mathematically  valid  only 
in  the  range  of  heads  for  which  the  truncated  Maclaurin  series  approxi¬ 
mates  the  exponential.  The  range  for  which  it  will  be  valid  depends 
upon  the  range  of  heads  over  which  the  empirical  coefficient  was  deter¬ 
mined.  Holler  calibrated  the  DRM  for  low-head  structures  and  conse¬ 
quently  the  model  yields  less  accurate  predictions  for  heads  greater 
than  20  ft. 

*  Unpublished  data,  C.  H.  Tate,  Jr.,  August  1978,  U.  S.  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  Mississippi. 
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Figure  6.  DRM-predicted  D.O.  versus  observed  D.O. 


Enid  Outlet  Works 


23.  Although  the  DRM  calibrated  for  low-head  structures  is  not 
rigorously  valid  at  high  heads,  large  predictive  errors  will  not  neces¬ 
sarily  result.  As  indicated  in  paragraph  16,  the  standard  error  of  pre¬ 
diction  with  the  DRM  (calibrated  for  low-head  structures)  and  the  EDM 
(calibrated  with  Enid  data)  was  0.9  and  0.7  mg/£,  respectively.  This 
occurs  because  the  final  D.O.  concentration  asymptotically  approaches 
saturation  as  the  head  is  increased.  Consequently,  at  large  heads,  the 
deficit  ratios  became  small,  which  results  in  relative  insensitivity  of 
the  final  D.O.  to  the  deficit  ratio.  The  absolute  error  in  D.O.  predic¬ 
tion  will  be  small;  however,  the  DRM  will  tend  to  predict  a  lower  final 
D.O.  than  actually  exists.  This  result  is  clearly  indicated  by  the 
plots  in  Figure  7  and  is  verified  by  the  comparison  of  predicted  and 
observed  results  in  Figure  5. 

24.  It  is  possible  to  calibrate  the  DRM  for  large-head  structures; 
however,  large  changes  in  the  experimentally  determined  coefficient  re¬ 
sult.  Calibrating  the  DRM  with  Enid  data  resulted  in  a  8"  of  0.200 
per  ft  as  compared  to  0.065  per  ft  for  low-head  structures.  As  indi¬ 
cated  in  Figure  8,  with  this  modification  the  DRM  and  EDM  predict 
essentially  equivalent  results  for  heads  of  50-70  ft.  This  is  also 
reflected  in  Figures  4  and  6,  which  compare  predicted  and  measured  D.O. 
for  the  EDM  and  DRM,  respectively.  DRM,  calibrated  for  large  heads, 
should  not  be  used  for  low-head  structures.  As  indicated  in  Figure  7, 
very  large  errors  would  result. 

Recommendations 

25.  The  EDM  should  be  used  instead  of  the  DRM  to  predict  D.O. 

The  EDM  with  an  escape  coefficient  of  0.045  per  ft  predicts  the  deficit 
ratio  for  low-head  structures  as  well  as  the  DRM  does.  Additionally 
the  EDM  predicts  D.O.  for  high  heads  (>25  ft)  without  modifying  the 
experimentally  determined  coefficient.  In  comparison,  large  changes  in 
the  experimentally  determined  coefficient  are  required  for  the  DRM  to 
achieve  equivalent  accuracy.  The  assumptions  upon  which  the  DRM  is 


26.  Additional  research  should  be  conducted  to  evaluate  the 
effects  of  flow  rate  on  gas  transfer.  Detailed  studies  of  the  hydraulics, 
turbulence,  and  energy  dissipation  related  to  gas  transfer  must  be  per¬ 
formed  to  clearly  identify  the  cause  and  effect  of  structure  design  and 
gas  exchange.  Projects  with  outlet  works  similar  to  the  structures 
examined  but  with  other  hydraulic  conditions  such  as  submerged  outlet 
portals  and  full  conduit  flow  should  be  examined  to  evaluate  gas- 
transfer  characteristics.  Since  reservoir  outlet  works  are  only  a 
portion  of  the  structures  used  to  release  water  from  man-made  lakes, 
other  types  of  structures,  such  as  spillways,  sluices,  and  hydropower 
structures  (Figures  9  and  10),  should  also  be  investigated  to  evaluate 
appropriate  models  with  these  gas-transfer  characteristics. 
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Predicted  by  SELECT  to  be  entering  outlet  works  prior  to  reaeration. 
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Appendix  B:  Example  Case  for  Prediction  of  Reaeration 
for  Gated-Conduit  Outlet  Works 


Synopsis  of  Procedure 

1.  Based  upon  the  analysis  presented  earlier,  the  D.O.  released 
from  gated  conduits  may  be  computed  if  the  following  conditions  are 
satisf ied : 

a_.  Flow  issued  from  under  the  bottom  edge  of  the  control  gate 

b^.  The  flow  rates  were  such  that  free-surface  flow  existed  in 
the  conduit;  i.e.,  the  conduit  was  not  flowing  full. 

c^.  The  outlet  portal  of  the  conduit  was  not  submerged  by  the 
tailwater . 

2.  The  data  required  are  the  following: 

a..  Lake  temperature  and  D.O.  profiles.* 

Forebay  and  tailwater  elevations. 

Structure  details,  i.e.,  outlet  configuration  and  location 
cL  Discharge  rate. 

_e.  Upstream  morphometry. 


3.  The  recommended  calculation  sequence  is  as  follows: 

a_.  Predict  the  flow  weighted  average  temperature  and  D  0. 
entering  the  outlet  works  with  the  computer  code  SELECT 
(Bohan  and  Grace  1973). 

b .  Adjust  the  escape  coefficient  for  the  temperature  with  the 
equation 


CT  ~ 


c20  i'022 


T-20 


where 

=  temperature-adjusted  escape  coefficient,  ft 

T  =  Predicted  release  temperature,  °C 

c?n  =  Escape  coefficient  at  20°C  (or  c9n  =  0.043 
U  per  foot)  ZU 


*  Either  observed  or  predicted  with  a  numerical  model  such  as  WESTEX 
(Loftis,  in  press).  (References  mentioned  in  this  Appendix  are  more 
fully  identified  in  the  References  section  at  the  end  of  the  main 
text. ) 
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Compute  the  oxygen  deficit  D  (mg/ Ji)  entering  the  outlet 
works  by 


D  =  C 


sat 


"SEL 


where 

oxygen  saturation  concentration  for  predicted 
release  temperature,  mg/ J l 

D.O.  entering  outlet  works  prior  to  reaera¬ 
tion,  predicted  by  SELECT,  mg/ Z 

d.  Calculate  the  difference  in  forebay  and  tailwater  eleva¬ 
tions,  Ah  ,  ft. 

e.  Predict  the  final  or  release  deficit  from 


Df  =  Di  exp  (-cT  Ah) 


where 

=  downstream  deficit,  mg/ i 

f.  Calculate  the  released  D.O.  concentration  from 


"rel 


-  C 


sat 


-  Dr 


Example  Case 


An  example  case  is  presented  below: 

a_.  The  temperature  and  D.O.  profiles  and  hydraulic  data  for 
Mystery  Lake  are  presented  in  Tables  B1  and  B2. 

_b.  These  data  are  input  to  SELECT  (presented  in  Appendix  C) 
and  the  following  predictions  are  obtained: 


T  =  18.34  °C 
CSEL  =  3-95  mg/ i 


_c.  The  temperature-corrected  escape  coefficient  is 
cT  =  0.045  (1.022) (T-20) 
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c  =  0.045  (1.022)(18,34  20) 

10.34 


=  0.043  ft 


and 

Ah  *  400  ~  340  -  60  ft 

d_.  The  upstream  deficit  is 

D.  -  C  -  3.95 
l  sat 

-  9.5  -  3.95 
=  5. 55  mg/£ 

e.  The  downstream  deficit  after  reaeration  is- 


Dj.  =  exp  (-cAh) 


=  5.55  exp  Q-0.043  (60)] 


=  0.4  mg/ £ 

f.  Released  D.O.  after  reaeration  is 


C  -|  =  C  -  D,. 
rel  sat  f 


=  9.5  -  0.4 
=  9.1  mg/ £ 
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Table  Bl 

Mystery  Lake  Profiles 


Elevation 

Depth 

Temperature 

D.0. 

ft 

ft 

°C 

mg/2 

400 

0 

28.0 

8.0 

395 

5 

27.6 

8.0 

390 

10 

27.3 

8.0 

385 

15 

26.0 

8.0 

380 

20 

25.0 

7.3 

375 

25 

21.0 

4.3 

370 

30 

14.0 

3.1 

365 

35 

12.0 

2.8 

360 

40 

10.1 

2.6 

355 

45 

10.0 

2.5 

350 

50 

9.9 

2.3 

Table  B2 

Hydraulic  and  Geometric  Data 


Discharge:  225  cfs 

Port  Location:  Elevation  372  ft 
2 

Port  Area:  16  ft 
Upstream  Morphometry: 


Elevation, 

ft 

Width,  ft 

400 

1500 

390 

1200 

380 

1900 

370 

800 

360 

600 

350 

100 

Forebay  Elevation: 

o 

o 

ft 

Tailwater  Elevation:  340  ft 
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Appendix  C 


SELECT  Example 


Cl 


1000 

J<J<J<MYSTERY  LAKE 

WITHDRAWAL 

1010 

J<J<J<FILES 

05 

06 

1020 

J<J<J<DATA  SETS 

1 

1030 

J<J<J<PRINT 

INPUT 

1040 

J<J<J<MYSTERY  LAKE 

EXAMPLE 

CASE 

1050 

J<J<J<ENGLISH 

1060 

J<J<J<TABLES 

1 

1070 

J<J<J<THICKNESS 

2.0 

1  OS  0 

J<J<J<INTERVAL 

2 

1090 

J<J<J<SURFACE 

400  . 

1100 

J<J<J<BOTTOM 

350  . 

1110 

J<J<J<NUMBER 

OF  WIDTHS 

6 

1120 

J<J<J<ELEVATION  WIDTH 

1130 

XX  J< 

400  . 

1500  . 

1140 

XXX 

390  . 

1200  . 

1150 

XXX 

380  . 

1  0  0  0  . 

1160 

XXX 

370  . 

800  . 

1170 

XXX 

360  . 

600  . 

1  IS  0 

XXX 

350  . 

200  . 

1190 

J<J<J<PORTS 

1 

1200 

xxxaREA 

16  . 

1210 

J<J<J<ELEVATION 

372. 

1220 

J<J<J<FLOW 

225. 
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